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We study the radiation from a collision of black holes with equal and opposite linear momenta. 
Results are presented from a full numerical relativity treatment and are compared with the results 
from a "close-slow" approximation. The agreement is remarkable, and suggests several insights 
about the generation of gravitational radiation in black hole collisions. 



I. INTRODUCTION 

The collision of two black holes is now being studied extensively via the techniques of numerical relativity [0. 
Collisions are of great importance as the most interesting source of gravitational waves that might be observable with 
interferometric detectors Q . The study is also of great inherent interest to relativity theory in that supercomputers 
allow us to investigate strong field gravity effects without symmetries which might preclude interesting or crucial 
phenomena. In dealing with such a daunting problem, useful checks, guidelines, and insights have been provided by 
analytical approximations, in particular by the close-limit approximation H. In principle, this method applies when 
the holes are initially very close together. In this case, the horizon is initially only slightly nonspherical and the 
spacetime that evolves outside the horizon can be treated as a perturbed single black hole. The highly nonspherical 
nature of the spacetime inside the horizon is causally disconnected from the exterior, and from the generation of 
outgoing gravitational waves. The exterior spacetime can be evolved forward in time from the initial data hypersurface 
with the linearized equations of perturbation theory. 

This method turns out to be remarkably successful The details of this success may give insights into the 

nature of collisions of holes. For holes that are initially momentarily stationary, the close-limit predictions of radiated 
energy and waveforms are quite good (i.e., in agreement with the results of numerical relativity) even when the initial 
horizon is highly distorted, violating the assumptions underlying the method. The close limit has been used by 
Abrahams and Cook B for the head-on collision of holes with initial momenta towards each other. This momentum 
causes horizons to form when the holes are at larger separation and makes the exterior spacetime more spherical, so 
it is not surprising that the close limit should be successful for these cases. Puzzling results emerge, however, when 
close-limit calculations are combined with Newtonian trajectories to estimate the radiated energy for initially large 
separations of initially stationary holes. The success of these estimates suggests, among other things, that to a large 
extent the role of the early weak-field phase of the evolution is to only to determine what the momentum of the holes 
will be when they start to interact nonlinearly. 

With that suggestion as one of our motivations, we consider here equal mass holes which are initially moving 
towards each other with equal and opposite momentum P. We analyze the problem with an approximation simple 
enough to allow insight, and we present, for comparison, the results of full numerical relativity for the same initial 
black hole configuration. In a certain sense, this study complements that of Ref. fij. The initial data sets being 
studied are representations of the same physical system; in Ref. j?J the data were "exact" (up to numerical error) 
solutions to the initial value problem, however, in the current study we have more control over the approximations 
implicit in the perturbative analysis. In Sec. II we present the general formalism for the problem and briefly discuss 
the full numerical solution. In Sec. Ill we describe an approximation based on the close limit and on slow initial 
motion. Results of both methods are presented and discussed in section IV. Throughout the paper we use units in 
which c = G = 1, and M represents the total ADM mass on the initial hypersurface. 

II. INITIALLY MOVING HOLES 

The initial value equations for general relativity are j|] , 



1 



V a {K ab -g ab K) = Q (1) 
3 R - K ab K ab + K 2 = (2) 

where g ab is the spatial metric, K ab is the extrinsic curvature and 3 R is the scalar curvature of the three metric. One 
proposes a three metric that is conformally flat g ab = ^gab, with g ab a flat metric, and <p A the conformal factor, and 
one uses a decomposition of the extrinsic curvature K ab = 4>~ 2 K ab . The constraints become, 

V a K ab = (3) 

V 2 0= -h- 7 K ab K ab , (4) 
o 

where V is a flat-space covariant derivative. 

In describing how (j^) and (Q) and the 3+1 evolution equations are solved numerically, it is useful to have at hand 
three different coordinate systems. Of greatest relevance to the numerical method are the Cadez coordinates, a 
system which is particularly well-suited for the collision of two black holes and which has been used extensively in 
numerical studies |],[l(| . These coordinates are spherical near the throats of both holes and in the asymptotic wave 
zone, so they simplify the application of both inner and outer boundary conditions. It is useful also to refer to two 
coordinatizations of the flat conformal three space: cylindrical coordinates p,z,(p, and the bispherical-like Misner |Tl| ] 
coordinates p,,r], ip. The fact that the problem is axisymmctric, of course, reduces the spatial computational grid to 
a two dimensional one. By choosing to consider only equal mass holes with equal and opposite momenta, we have a 
further symmetry which reduces the size of the computational grid to a quadrant, (jp = 0, z > 0). We characterize 
the separation of the holes with the Misner parameter fj,Q, and construct the coordinate grid independently for each 
choice of Details of the grid computation are given in Refs. fi~^ , |i~3l| . 

To solve the momentum constraint (||) we follow the prescription of York and coworkers |l4|] and Cook |l5|,[l(| . This 
starts with a solution to (||) that represents the momentum of one hole, 

KIT = [2P (a n b) - (S ab - n a n b )P c n c ] . (5) 

Here the hole is associated with some point in the flat conformal space, r is the vector from that point, and n is the 
unit vector in the r direction. The next step is to modify (|J) to represent holes centered at z = ±coth/i , the centers 
of the circles /i = ±/io in the conformally flat metric. Since the momentum constraint (|^) is linear, one can simply 
add two expressions of the form (||): 

K\J° = K°f c (*-»*- coth fi ) + K°f c {z^z + coth ^o,P^-P) ■ (6) 

For convenience, the initial data is forced to obey an isometry condition, i.e., we operate on the momentum constraint 
solution with a reflection procedure equivalent to adding image charges in electrostatics. The result of this procedure 
is to create a solution which corresponds to two identical asymptotically fiat universes connected by two Einstein- 
Rosen bridges. The nature of this symmetrization process, and the boundary condition it provides for (Q), affects the 
mass of the holes being represented. Cook Jl6|,[l7j has also used this approach to develop codes to compute symmetric 
initial data solutions for axisymmetric and full 3D data. 

The Hamiltonian constraint is solved by linearizing equation (Q) around a solution (j>\ so that (f> — <f>i+6<p, discretizing 
to second order the resulting linear elliptic equation, solving the matrix equation for 5<f) with a multigrid method, 
then iterating the procedure until a convergence tolerance of §<fi/<f>i < 10 is achieved. It has been verified that for 



K a b — 0, the solution for <j> converges quadratically with cell size to the time-symmetric Misner data |11|. 

The symmetrized initial data for <fi and for K ab are now used as the starting point for numerical integration. The 
evolution employs maximal time slices and the shift is determined by an elliptic condition that forces the 3-metric (in 
Cadez coordinates) into diagonal form pj . The numerical errors inherent in the method (to be described elsewhere) are 
similar to those in Ref. 0] . We have verified that the convergence rate for the total radiated energy scales quadratically 
with grid spacing and that differences in the dominant £ = 2 waveforms are on the order of a few percent at the grid 
resolutions used here. The errors are small on the scale of Fig. ||, and do not affect any conclusions to be drawn 
from that figure. The methods used for the numerical evolution are described in detail in Ref. [^0|; we modified only 
slightly the code described there for evolving the time symmetric Misner data. 



III. APPROXIMATION METHOD 



The close-limit approach can be applied to the Cook [|T7| initial data, as has been done in Ref. Q. But the Cook 
initial solution is numerical. To facilitate insights we make a further approximation. We assume that the black holes 
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are initially close, and that the initial momentum P is small. Our solution for the extrinsic curvature K a f, is 
(^), the simple superposition (without symmetrization; this effect will be discussed later) of two one- hole solutions. 
We denote by n + and ri~ the normal vectors corresponding, respectively, to the one hole solutions at z = +L/2 and at 
—L/2, and we define R to be the distance to a field point, in the flat conformal space, from the point midway between 
the holes. For large R, the normal vectors n + and ft" almost cancel O]. More specifically n + = —-nT + 0(L/R). A 
consequence of this is that the total initial K ah is first order in L/R, and its (R,9,ip coordinate basis) components 
can be written as 
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In addition to being first order in L, the solution for K ab is first-order in P and therefore the source term on the 
right in the hamiltonian constraint (f|) is quadratic in P. If we limit ourselves to a solution to first order in P we can 
ignore this quadratic source term. (In Sec. IV, a more thorough discussion will be given for this step of ignoring the 
source term.) Without the source term the hamiltonian constraint reduces to the zero momentum case, the Laplace 
equation. The symmetric solution to this (i.e., the solution for two identical asymptotically flat universes) is the Misner 
solution jn]], and this is the solution we take. The Misner geometry is characterized by a dimensionless parameter 
Ho which describes the separation of the throats. We must, of course, choose /io appropriate to the parameters of the 
extrinsic curvature we are using. We choose therefore a Misner geometry characterized by the same value of L as in 
(0). Since L there represents not the physical distance, in any sense, between the holes, but the formal distance in 
the conformally flat space, we choose a Misner geometry with the same value of L in the conformally flat part of the 
Misner metric. The relationship of L to /io is (see, e.g., |Q) 
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This completes the description of the initial data to first order in L and to first order in P (the close, slow approxima- 
tion). We now view the spacetime exterior to the horizon as a perturbation of a single Schwarzschild hole described 
in standard Schwarzschild coordinates t, r, 0, <p. Even-parity perturbations are then described by a Zerilli function if). 
According to the general prescription given in reference || the value of ip on a t = initial hypersurface is found from 
the initial value of the three geometry. Our initial geometry, to first order in P, is exactly the same as the zero P 
solution in reference [Q, where the Zerilli function is denoted V'pert, and is given in eq. (4.29), along with (4. 10), (4. 27) 
and (4.28). In that reference it is shown that in the close limit, the quadrupole contribution dominates, with contri- 
butions for I > 2 higher order in the separation parameter. Here we shall consider only the 1 = 2 contribution, and 
shall denote the Zerilli function, corresponding to this Misner (i.e., P = 0)problem, as tpMis(r, t). 

The initial value of ip, the time derivative of the Zerilli function, follows from the extrinsic curvature as explained 
in The extrinsic curvature is given, in our approximation, by multiplying K ah in (0) by the squared reciprocal of 
the conformal factor for the Schwarzschild geometry, 0schw = 1 + M/2R. We must map the coordinates of the initial 
value solution to the coordinates for the Schwarzschild background. To do this, we use the same mapping used for 
the initial value of if> in Q: we interpret the R of (fj]) as the isotropic radial coordinate of a Schwarzschild spacetime, 
and we relate it to the usual Schwarzschild radial coordinate r by R = (y/r + \Jr — 2M) 2 /4. From this we arrive at 
the following expression for the (Schwarzschild coordinate basis) components of the extrinsic curvature: 
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Here we have used the fact that 
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From (|^), which contains both monopole and quadrupole parts, we can project out the 1 = 2 part and read off the 
initial value of the time derivative of the Zerilli function to be 



2APLJ1 - 2M/r [4x ' 3M, 
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(11) 
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FIG. 1. The anticorrelation between the two components of the perturbation, i/>Mis and ipMom, leads to the negative coefficient 
in the energy vs. momentum relation. 

Along with ip\t=o — ipMis(r,t = 0), this completes the specification of the Cauchy data for ip. 
Given this Cauchy data, the time evolution is obtained by evolving the Zcrilli equation, 

d 2 ^/dt 2 - d 2 ip/drl + V{r)ip = , (12) 

where = r + log(r/2Af — 1) and the Zcrilli function ip is a coordinate invariant combination of the perturbed metric 
coefficients; the 1 = 2 "potential" V(r) can be seen in reference 
The evolved ip can be decomposed into two components 

Ip = IpMis + IpMom ■ (13) 

The first term is the solution of ( |l2|) for cauchy data ip = ipMis(r, t = 0) and ip = at t = 0. The second term is the 
solution for t = cauchy data ip = 0, and with ip given by (pd|). The two contributions are respectively zero order in 
P and first order in P; the decomposition then represents a separation into parts of ip due to the masses, and to the 
momenta. The radiated energy is given by [^) 

1 f°° . 

E = — — / ip 2 dt , (14) 
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and can be written, in terms of the decomposition above, as: 

E = yj 1pM iB dt + 2 J IpUistpMom dt + J ^Mom dt J ■ ( 15 ) 

The first term gives the same result as in the momentarily stationary case; it is simply the radiation for the Misner 
initial geometry, as computed in reference . The second term is linear in the momentum of each hole. The coefficient 
of it is given by the "correlation" of ^Mis and V'Mom- As can be seen in Pig. this "correlation" integral is negative. 
The anticorrelation is compatible with previous simulations done by Ref. |7(] using numerical initial data (see figures 
3a,b in their paper). This means that for small values of P, the radiated energy decreases with increasing momentum. 
The effect is clearly visible in Fig. |^ where we show the radiated energy as a function of the momentum, 

Note that the first term is simply a function of the Misner parameter /io- The second term depends on /iq, but 
also depends on L and P. We can write L in terms of /io with to express all dependencies in (|l^) only in terms 
of /ito and P/M. With the correct numerical factors we get the final result of the close-slow approximation, a simple 
formula for the radiated energy simply and explicitly expressed in terms of the parameters of the collision: 
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FIG. 2. Energy as a function of initial momentum. Here E is the energy radiated during coalescence, P is the initial 
momentum, and M is the initial ADM mass. Curves shown are for fixed parameter /io , corresponding to separation of the holes 
in conformal space. The curves clearly show the "dip" effect, and the good agreement, even for large values of the momentum. 



In Fig. |2|, we plot the radiated energy computed from (Jlfj) for several values of initial separation /io, and for a wide 
range of P/M . On this plot, also, are presented the results for radiated energy from numerical results computations 
which make no approximations. The agreement between the numerical results and the results of the approximation 
is remarkably good, even at rather large values of P/M. 



Two features of Fig. g stand out. The first is "momentum dominance": the radiated energy is dominated by the 
third integral in ( |l5| ) unless the momentum is very small. The second obvious feature is that the approximation 
method works very well even for sizeable values of P/M. 

To understand the implications of these features, let us start by reviewing the difference between the exact, nonlinear 
numerical computation, and the approximation scheme of Sec. III. In the exact method we start with an exact solution 
to the initial value equations described by two parameters, one a dimensionless measure of the separation of the holes, 
the other a dimensionless measure of the momentum. The process of generating the solution consists of four steps: (i) 
One starts with a very simple prescription for K a t, constructed by superposing two solutions of form (||) corresponding 
to two coordinate positions in the conformally flat space, (ii) Equation (Q) is then solved for the conformal factor 
and hence for the three geometry, (iii) The solution for the extrinsic curvature and the initial geometry is then 
"symmetrized" by an iterative process equivalent to adding image charges, (iv) This solution is numerically evolved 
off the initial hypersurface with the full nonlinear Einstein equations. 

By contrast, the steps for the approximate solution are: (i) The (conformal) extrinsic curvature is taken to be the 
unsymmetrized superposition of two contributions with the form of (|^). (ii) The conformal factor, and therefore the 
three geometry, is taken to be the symmetrized solution corresponding to throats located at the same points in the 
conformally flat space as the points in K a t. (iii) This approximate initial data is then treated as initial data for the 
nonspherical perturbations of a Schwarzschild hole, and the perturbations are evolved with the linearized Einstein 
equations. 

The difference in evolution off the initial hypersurface (full Einstein equations in one case, linearized equations 
in the other) is not a major source of error in the interesting cases, those with high momentum. As momentum 



where «2, as defined in Ref. ||, i 



is 




(17) 



IV. RESULTS 
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increases, the location of the horizon in the initial geometry moves outward. The high momentum cases, therefore, 
correspond to throats which, on the initial hypersurface, are well inside an all-encompassing horizon. This is the 
situation in which the "close-limit" approximation method should work very well. It is also not surprising that no 
large error is introduced by the failure, in the approximation method, to symmetrize the extrinsic curvature. One 
way of understanding this is to note that ipyi om lacks the "image" contributions needed for symmctrization. These 
images only influence the form of K a b very close to the holes. As the separation between the holes gets smaller the 
horizon moves further from the throats and the effect of the images on K a b outside the horizon diminishes. We have 
checked numerically that the difference between the symmetrized and unsymmetrized K a b, for all cases considered, is 
negligibly small outside the horizon. 

These two aspects of the approximation method rely on the throats being "close" in some sense, an approximation 
that seems well justified. What remains to be explained is how the slow-limit approximation does such a good job 
of approximating the very "unslow" correct initial data. We must also justify the apparent inconsistency in how the 
approximation scheme deals with orders of P. In the computation of ip the scheme explicitly omits corrections of order 
P 2 in (|4|). Formally, then, we should only be able to keep terms of first order in P in (JlJ). But it is the apparently 
inconsistent P 2 terms, of course, which dominate at most points in Fig. || ("momentum dominance" in generation of 
radiation). Not only do the P 2 terms agree with the results of numerical relativity, but the agreement remains good 
for rather high values of P/M . This raises the question: just what momentum contributions has our approximation 
really omitted? 

The momentum enters into the construction of the initial data in only two direct ways. First, it is an overall scaling 
parameter for K a b- The expression in (^) is an approximation for small L, but it is exact in P. The process of 
symmetrizing does not change this. Up to a conformal factor, then, the extrinsic curvature is exactly linear in P. 
Second, P enters the determination of the conformal factor through (^). The success of the slow approximation must 
be directly ascribed to the relatively unimportant role played by the right hand side of (Q). 

Further work will be needed for a real understanding of this, but some reasonable speculations can already be made. 
Due to momentum dominance the details of the initial three geometry are not crucial, so any quadrupolar distortion 
induced by K a b at large P will be insignificant compared to the radiation generated by the extrinsic curvature. The 
"slow" approximation, of course, is not perfect; at sufficiently high momentum it begins to fail. We speculate that 
the reason for this failure is not primarily due to K a b generating quadrupolar distortions of the initial three geometry. 
Rather, it is the effect of that source on the monopole part of the conformal factor, and hence on the ADM mass "M," 
that is used to scale physical quantities. When we do a comparison in Fig. || between the numerical relativity results 
and those of the approximation, we are comparing two cases for the same [Iq (i.e., the same coordinate separation in 
conformal space) and for the two cases we compare E/M at a given value of P/M. We are therefore placing on an 
equal footing the true value of M in the numerical relativity solution, and the P — > value of M in the approximation. 
It should be possible, in principle, to correct for this and, in effect, reduce the approximation to one in which we have 
only ignored the quadrupolar part of the source in (Q). 

The present results greatly help us to understand the success of the results of Ref. §. That success seems to 
require two things about the generation of gravitational radiation in collisions from large distances: (i) There must 
be negligible radiation during the early motion, when the holes are in each other's weak field region, (ii) The only 
important consequence of the early, weak-field, motion must be to give the holes momentum when each reaches 
the strong field region of the other. The first requirement is relatively easy to check. In Fig. || we plot radiated 
energy, computed by methods of numerical relativity, as a function of time, first for initial data representing two 
black holes falling from large separation. (The oscillations are due to the fact that almost all the energy comes off 
as "quasinormal ringing" of the final hole formed.) We also show the result of a second calculation. Cook [[l7j initial 
data are taken corresponding to the separation and momentum that the black holes would have after falling to a fairly 
close separation. A comparison of the curves verifies that the early stage of motion does not produce a significant 
contribution to the total outgoing radiation. 

Our present results, and in particular momentum dominance, strongly support the second requirement for the 
success of the ideas of Ref. 0. Since ipMom is the source of essentially all the radiation, one can see that what is 
important about the early stages of the coalescence is only the development of extrinsic curvature. This does not, of 
course, explain why there seems to be insensitivity to the details of the extrinsic curvature. (Surely, the Bowen-York 
extrinsic curvature, symmetrized or not, is not actually the extrinsic curvature that evolves from earlier stationary 
conditions. Yet, it seems to be adequate to give good predictions.) A more satisfactory answer to this question means 
that we must understand the relationship between data on an initial hypersurface and how this evolves to data on 
subsequent hypersurfaces. We must also understand the importance of confining ourselves to conformally flat data 
on hypersurfaces. Progress on these questions will probably require comparable results from four distinct classes of 
initial data sets. These are (a) Misner data with large hole separation, (b) the non conformally-flat data with close 
holes that evolves from (a), (c) boosted conformally-flat data with close holes, and (d) boosted conformally-flat data 
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FIG. 3. Radiated energy as a function of time for two different initial value sets. The first is for equal mass holes falling 
from rest at p = 2.2. The second shows the result of a boosted collision starting from a separation parameter /xo = 1.406 and 
a momentum P/M = 0.23. The second set of initial data can be considered to be an approximation to a constant t slice of the 
spacetime that evolves from the first set. The time scale is the (flat space) retarded time with zero corresponding to the time 
of apparent horizon formation. The energy scale of the Ho = 2.2 curve has been offset to zero at the time of apparent horizon 
formation. 

in the close-slow approximation. In addition, one requires reasonable measures of physical separation and momentum 
so that correspondence can be drawn between disparate initial data sets. 

There is strong motivation for carrying out such studies. The results so far achieved, both by numerical relativity 
and with the close and the slow approximation, are limited to head-on collisions. The situation of astrophysical 
interest, of course, is very different: the coalescence of orbiting holes. If the last few orbits in a coalescence are to 
be studied with numerical relativity, it will be crucial to understand what initial data are to be used to start the 
computation. Studies with the head-on collision provide a useful starting point to understanding the sensitivity of 
the radiation generation to the details of the initial data. 

A rather different, and more speculative, motivation for a better understanding of these issues, is the hope that 
our approximation methods might be as successful with orbital problems as with head-on coalescence. These results 
might provide "easy" approximate answers over a reasonable range of orbital coalescences, and may therefore serve 
as a guide to the numerical studies. 
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